Semiclassical theory of weak antilocalization and spin relaxation 

in ballistic quantum dots 
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We develop a semiclassical theory for spin-dependent quantum transport in ballistic quantum dots. 
The theory is based on the semiclassical Landauer formula, that we generalize to include spin-orbit 
and Zeeman interaction. Within this approach, the orbital degrees of freedom are treated semiclas- 
sically, while the spin dynamics is computed quantum mechanically. Employing this method, we 
calculate the quantum correction to the conductance in quantum dots with Rashba and Dresselhaus 
spin-orbit interaction. We find a strong sensitivity of the quantum correction to the underlying 
classical dynamics of the system. In particular, a suppression of weak antilocalization in integrable 
systems is observed. These results are attributed to the qualitatively different types of spin relax- 
ation in integrable and chaotic quantum cavities. 

PACS numbers: 03.65.Sq, 71.7D.Ej, 73.23.-b 



I. INTRODUCTION 

Guided by the vision to incorporate spin physics into 
the far-advanced semiconductor (hetero)structure tech- 
nology, semiconductor-based spin electronics (see, e.g., 
Ref. [lj) has developed into a prominent branch of present 
spintronics research. In this context spin-orbit (SO) in- 
teractions have recently received considerable attention 
since they give rise to rich spin dynamics and a vari- 
ety of spin phenomena in nonmagnetic semiconductors. 
Though SO couplings have been a subject of continu- 
ous research throughout the last decade a 2 i 3 i 4 i 5 i 6 i 7 i 8 , there 
is presently a revival in investigating SO effects owing to 
their role in spin transistors&Afi, spin interferometersii*i£, 
spin filtersi&i^, and spin pumps 1 ^^, to name only a few 
examples. Furthermore, most recently the intrinsic spin 
Hall effect in a SO-coupled system has caused an in- 
tense and controversial discussion in the literature^. Fi- 
nally, in spin-based quantum computation SO-induced 
spin relaxation effects may play a rolei^. 

The interplay between spin dynamics and confine- 
ment effects is particularly intriguing in quantum trans- 
port through low-dimensional devices at low tempera- 
tures where quantum coherence effects additionally arise. 
There exist two prominent experimental probes for SO 
effects in quantum transport: (i) characteristic beat- 
ing patterns in Shubnikov-de Haas oscillations in two- 
dimensional electron gases with tunable SO coupling, 
controlled via a back-gate voltage 20 - 2122 - 23 , and (ii) 
weak antilocalization^i 2 ^ 2 ^ (WA), an enhancement of 
the magnetoconductance at zero magnetic field owing 
to spin-dependent quantum interference effects. Since 
systems without SO coupling exhibit weak localization 
(WL), i.e., a reduction in the magnetoconductance, the 
appearance of WA allows conclusions to be drawn on the 
SO strength. While WA is fairly well understood for dis- 



ordered bulk system 
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relaxation phenomena in confined systems. These mea- 
surements are focussed on the interesting inter-relation 
between quantum confined orbital motion and spin evo- 
lution and relaxation in clean ballistic quantum dots. In 
these systems, the elastic mean free path is exceeding 
their size, and impurity scattering is replaced by reflec- 
tions off the system boundaries. 

Corresponding efforts in treating SO effects on spec- 
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, spin relaxatio: 
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and the interplay be- 



in recent experiments using 



ballistic bismuth 3 ^ and GaAs^i cavities, WA has been 
employed to study SO-induced spin dynamics and spin 



tween SO and Zeeman coupling 3 ^ 4 ^ in quantum dots 
have also been made on the theoretical side. SO-induced 
WA in ballistic quantum dots has been studied using 
random-matrix theory 41 ^ (RMT) and semiclassical ap- 
proaches^ 3 ^. While RMT approaches are restricted to 
quantum dots with corresponding chaotic classical dy- 
namics, the semiclassical transport theory comprises a 
much broader class of systems, including integrable con- 
finement geometries. Related semiclassical techniques 
have also been applied to spin transmission 45 and spin 
relaxation 4 ^ in quantum dots. 

The purpose of the present article is a detailed ex- 
position and extension of the semiclassical methods of 
Refs. an d 13 The theory to be discussed here 
unifies two subject areas: the semiclassical descrip- 
tion of WIj 4 ^ 4 ^ 4 ^ and the semiclassical treatment of 
SO interaction 50 i 51 i 52 i 53 i 54 i 55 i 56 . Compared to the earlier 
works 43 44 - 46 on the subject, here we give special atten- 
tion to the differences in spin relaxation along open and 
closed trajectories, analyze the interplay between Rashba 
and Dresselhaus interaction, and report on the full quan- 
tum calculations of spin-dependent transmission and re- 
flection. 

This paper is organized as follows: In the introductory 
Sec. m using path integrals with spin coherent states 
we deduce a spin-dependent semiclassical propagator and 
the corresponding Green function. Our main analytical 
results are presented in Sec. IIIII There, on the basis 
of Green functions a semiclassical approximation to the 
Landauer formula with spin is derived. The semiclassi- 
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cal Landauer formula is then applied to chaotic quantum 
dots, whereby the quantum corrections to transmission 
and reflection are calculated. In Sec. |^ we discuss how 
WA is related to the spin evolution. We define mea- 
sures for spin relaxation and consider, as an example, 
the spin relaxation in diffusive systems. In the following 
two sections the general theory is applied to chaotic and 
integrable quantum dots with Rashba and Dresselhaus 
SO interaction. In Sec. E] a detailed numerical study of 
the spin relaxation is followed by an analysis of the limit 
of slow spin dynamics (i.e., extremely weak SO coupling). 
Additionally, we examine a gauge transformation of the 
spin-orbit Hamiltonian that can be carried out in this 
limit. The dependence of the quantum corrections to 
transmission and reflection on the SO-coupling strength 
and magnetic field (Aharonov-Bohm and Zeeman contri- 
butions) is presented in Sec. I VII There some of the semi- 
classical numerical results are compared with full-scale 
numerical quantum calculations. 



II. SPIN-ORBIT INTERACTION IN A 
SEMICLASSICAL THEORY 

In this preparatory section we construct a spin- 
dependent semiclassical propagator— and related Green 
function. It fully describes the system at a given level 
of approximation and, thus, can serve as a starting point 
for our derivation of a semiclassical Landauer formula for 
systems with SO and Zeeman interaction fSec. 1111(1 . 

In the spinless case, the semiclassical propagator is 
conventionally obtained from the path-integral represen- 
tation of the exact propagator—. After the stationary- 
phase evaluation, which is valid in the semiclassical limit, 
the classical trajectories are selected from all the paths 
in the integral. In order to include spin into the path in- 
tegral, a continuous basis of spin states is required. The 
spin coherent states represent such a basisS*— 

Following Ref. |5^, we define a coherent state of spin 



i,l,..-by 



|C> = (l + |C| 2 )- s exp(Cs + )| C 



(1) 



where ( is a complex number that labels the state, s+ 



is the spin operator, and 



a) are the eigenstates 
, , s. To each ( corre- 



of s z with eigenvalues a — — s, 
sponds a three-dimensional unit vector n(£) = (£|s|£)/s 
that denotes the spin direction. It is easy to show that 
£ is a stereographic projection from the unit sphere cen- 
tered at the origin onto the plane z — 0. The projection is 
given by {n Xl n y ,n z ) i — > (Re C, — ImC 0)j where, in par- 
ticular, the south pole is mapped to Q = 0. In general, 
coherent states have the minimal uncertainty of s among 
all spin states and are characterized by three real param- 
eters: the direction n and an overall phase. (Hence, any 
state of spin 1/2 is coherent.) In the current definition, 
the phase is assigned to each n by Eq. , but other phase 
assignments are possible. Note that the phase of the state 
|C = oo) oc \<r = s) with n = (0, 0, 1) is not well defined. 



However, this manifestation of the fundamental problem 
of phase assignment^ does not pose a difficulty in our 
case, since the final results will be transformed to the 
\o~) representation using the projection operators |C)(CI- 
The states are normalized to unity, but, obviously, 
not mutually orthogonal (no more than 2s + 1 states of 
spin s can be mutually orthogonal). Nevertheless, having 
the property of resolution of unity, 



|C)(C|d/i(0 = i, dfi(Q 



2s + 1 



tt(i + ICI 2 )' 



:d 2 (, (2) 



they form an (overcomplete) basis and enable a path- 
integral construction. 

Let us consider a rather general case of a system with 
Hamiltonian linear in the spin operator s: 



H = H°(q,p) + hs-C(q,p). 



(3) 



Here q and p are the d-dimensional coordinate and mo- 
mentum operators, respectively, H° (q, p) is the spin- 
independent Hamiltonian, and hs ■ C(q, p) describes the 
SO interaction and the Zeeman interaction with an exter- 
nal (generally inhomogeneous) magnetic field. Utilizing 
the idea of Refs. |53| and l55l we express the propagator 
in the combined coordinate and spin-coherent-state rep- 
resentation in terms of the path integral 

C/(q 2 ,C2,qi,Ci;T) = (q 2 ,C 2 |e~ (l/ri)AT |qi,Ci) 



The integration is performed over the paths 
[q(^)> P(t)i C(t)] m the spin-orbit phase space con- 
necting (qi,Pi,Ci) to (q2,P2,C 2 ) in time T with 
arbitrary pi and P2. The integration measures are 
defined by 

-D/4C] = [[ d^(aW), (5) 



(2irh) 



where tj — jT/n. The Hamilton principal function 
W = W° + hsW 1 consists of two contributions: the usual 
classical part, 



W°[q,p;T]= f dt[p-q-H (q,p) 
Jo 



(6) 



and the spin-related part, 



W 1 [q, P ,C;T] 



dt 



cc* - c*c 

^(i + ICI 2 ) 



n(C)-C(q,p) 



(7) 

Now we can separate the integration over the spin paths 
in Q, thereby representing the propagator 



c/(q 2 ,C 2 ,qi,Ci;7 1 ) 

V[q]V[p] 



(27rft)° 



tf [q ,p](C2,Ci;T) exp -W°[q,p;T] (8) 
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with 

^ [q , p] (C2,Ci;T)= /'x> Ai [C]ex P {*aW 1 [q,p > C;T]}. (9) 

Clearly, if[q lP ] (C2 , Ci > T) is a propagator of a system with 
the time-dependent Hamiltonian i?[ qp ] (t) = fi,s-C[ qp ](i), 
where Cr q>p ](t) = C(q(t), p(t))is calculated along the 
path [q(i), p(t)] of the integral jHJ. This Hamiltonian de- 
scribes a spin, precessing in the time-dependent magnetic 
field C^p]^)^ Expression for if[ q . p ] (C2 , Cl! T) can 
be integrated explicitly 59 , yielding the usual spin propa- 
gator in the basis of coherent states (Appendix El . 

We proceed by evaluating the path integral © in the 
semiclassical limit W° ^S> h. The integration simpli- 
fies considerably, if the spin-dependent Hamiltonian is 
treated as a perturbation, i.e., when 



fo|C(q,p)| « |ff°(q,p)|. 



(10) 



This condition, assumed for the rest of the paper, is 
usually fulfilled in experiments based on the semicon- 
ductor heterostructures. According to this requirement, 
the spin-precession length must be much larger than 
the Fermi wavelength, however it can be smaller, of or- 
der, or greater than the system size. The semiclassical 
and perturbative regimes can be implemented simulta- 
neously^i*^ 2 . by formally letting h — ► and keeping all 
other quantities fixed. Then the phase of the integrand 
in Eq. JHJ) is a rapidly varying functional, which justifies 
the use of the stationary-phase approximation. It is cru- 
cial that ^[ q , P ] (C2, Cii T) does not depend on %, i.e., it 
is a slowly varying functional, and, therefore, its effect 
on the stationary trajectories can be neglected. Thus, 
the stationary trajectories are the extremals solely of 
W°[q, P',T], which means that they are the classical or- 
bits of the spinless Hamiltonian H . The resulting semi- 
classical propagator, 



[/ sc (q 2 ,C 2 ,qi,Ci;T) 

= ^x 7 (C 2 ,Ci;r)c 7 exp{^w 7 (q 2 ,q i; r)j, (11) 



is a sum over all classical trajectories 7 = [q 7 (i), p 7 (i)] 
of time T from qi to q 2 . The prefactor C 7 , arising from 
the stationary-phase integration, is the same as in the 
spinless case££: 



exp (- 



det 



a/3 



9 2 W0(q 2 ,qi;T) 



1/2 



where u 1 is the Maslov index 



dq%dqf 

(12) 

Although the classical 
trajectories are not affected by the spin motion, the 
reverse is not true. Indeed, the spin propagator K 7 , 
computed along the classical trajectories, describes the 
spin evolution in the effective magnetic field C 7 (i) = 
C(q 7 (t), p 7 (t)) generated by these trajectories. 



The semiclassical Green function is given by the 
Laplace transform of U SC (T) to the energy domain E, 



/>OC 

G(E) = -- dTe^ E+i0+ ^ T / h U BC (T), 
h Jo 



(13) 



evaluated in the stationary-phase approximation. As be- 
fore, K 7 (T) does not modify the stationary-phase condi- 
tion, and the theory without spin can be applied. More- 
over, using the resolution of unity the spin propaga- 
tor can be transformed to the usual \a) basis. Finally, 
we obtain 

GW(q2,qi;-E) 

= ^(# 7 ) ffV ^ 7 exp{~S°(q 2) q 1 ;£)l (14) 



with a, a' — —s,...,s. In Eq. Ijl4|l . 7 is a classical 
trajectory of the Hamiltonian H° — E with the action 
5° = / p - dq and time T 7 (E) = dS^/dE. K 7 (t) is 
the operator form of the spin propagator [Appendix 
Eqs. (M} and iftfijl]. and [K^) (t , a = {a'\K^ (T 7 (E)) \a) 
is its matrix element. The prefactor is given by 



-if sgn(dT T /<lE) 



(15) 



and C 7 is expressed in terms of the derivatives of 
5°(qa,qi;£) (Ref.|H|. 



III. SEMICLASSICAL LANDAUER FORMULA 
WITH SPIN 

The semiclassical Landauer formula with spin, derived 
below, is the main analytical result of this paper. It 
forms the basis for the subsequent semiclassical treat- 
ment of the spin-dependent transport in two-dimensional 
systems. 



A. Derivation of the formula 

We start from the standard (quantum) Landauer for- 
mula, that relates the conductance (e 2 /h)T of a sam- 
ple with two ideal leads to its transmission coefficient T 
(Ref. l63|) . Assuming that the leads support N and N' 
open channels (not counting the spin degeneracy), re- 
spectively, the transmission can be expressed as the sum 



N' N s 

T = |incr',mcr| 
n—1 m— 1 er', a— — s 



(16) 



Here t ncr > ma is the transmission amplitude between the 
incoming channel \ma) (with spin projection a) and the 
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outgoing channel \na') belonging to different leads. We 
shall also consider the reflection coefficient 



K = 



N s 

E E 

n,m=l a' u= — 



' ncr' ma 



(17) 



where the reflection amplitude r n(J / ma is defined for the 
channels of the same lead. The transmission and reflec- 
tion satisfy the normalization condition 



T + K=(2s + 1)N 



(18) 



that follows from the unitarity of the scattering matrix. 

Consider, as a model for a (large) quantum dot, a two- 
dimensional cavity (billiard) with hard-wall leads. The 
particle in the cavity is subjected to the SO and Zee- 
man interaction of the form ©. Semiclassical expres- 
sions for the transition amplitudes in the spinless case 
were derived in Refs. and by projecting a semi- 
classical Green function onto the lead eigenstates, while 
integrating over the lead cross sections in the stationary- 
phase approximation. For a particle with spin we im- 
plement this procedure using the semiclassical Green 
function Q14p. In the semiclassical limit of large action, 
5° ^> h, the spin-propagator element (X 7 ) (T / (T does not 
shift the stationary point. In the resulting expression, 



"n<7 ,mcr 



E 

7(n,m) 



exp ( -5 7 



(19) 



the only effect of spin is to weight the contribution of each 
trajectory in the sum with the respective matrix element 
of K 1 . In Eq. iJTSj l j(n, m) is any classical trajectory of 
energy E that enters (exits) the cavity at a certain angle 
©to (@fi) measured from the normal at a lead's cross 
section^ The angles are determined by the transverse 
momentum in the leads: sinG^ = rhir/kw and sin0 s = 
mr/kw', where k is the wavenumber and w and w' are 
the widths of the entrance and exit leads. The action for 
a trajectory of length L 7 is 5 7 = hkL-y. The prefactor is 
given by 



sgn(n) sgn(m) 



2 ^w' jcoseRCOse^M^I 172 



x exp 



ik (sin Qfh V - sin fi y') 



7T 

1 2 1 ih 



(20) 



where is an element of the stability matrix (as de- 



fined, e.g., in Ref. |6j), y (y') is the coordinate on the 
lead's cross section at which the orbit 7 enters (exits) the 
cavity, and (tx 7 is the modified Morse index 48 . Substitut- 
ing the sum (|19|l and the corresponding result for r na i iTna 
in Eqs. (|16fl and i|17fl we obtain the semiclassical approx- 



imation for the total transmission and reflection: 



(T,TZ) = 

E E E 

nm «-y(n,m) 7 / (n,m) 



M 7l7 '^4 7 ^4 7 / exp 



(S 1 — Sy ) 



(21) 



Here in the case of transmission (reflection) the paths 7 
and 7' connect different leads (return to the same lead). 
In this expression each orbital contribution is weighted 
with the spin modulation factor 



Tr(X 7 i^,) 



(22) 



where the trace is taken in spin space. Equations (|21|) and 
(|22|l generalize the semiclassical Landauer formulali 4 ^ to 
the case of spin-dependent transport. 



B. Leading semiclassical contributions for a 
spinless particle 

In the semiclassical limit the phases in Eq. (|21|l are 
rapidly varying functions of energy, unless 7 and 7' have 
equal or nearly equal actions. Therefore, if one calculates 
the transmission and reflection averaged over a small en- 
ergy window, most of the terms in the double sum will 
vanish. In the following, we review the leading contribu- 
tions for a spinless system (A4 7l7 ' = 1) with time-reversal 
symmetry: 

(i) The classical part consists of the terms with 7' = 7 
(Ref. l65j) . Their fast- varying phases cancel (including 
the phase in the prefactor). For a classically ergodic (in 
particular, chaotic) system one finda^ 



r (o) _ 



NN' 



N' 



n 



(0) 



N 2 



N + N' 



(23) 



(the superscript refers to zero spin and zero magnetic 
field). This result can be obtained using the sum rule^S 

]T lAyfSiL-L^^iN + N'y'PLiL). (24) 

7(n,m) 

It implies that the length L of the classical trajectories 
is distributed according to 



P L (i)^-^exp f-A 



(25) 



in other words, the probability for a particle to stay in 
an open chaotic cavity decreases exponentially with time. 
The average escape length is given by 



= ttA c = fcA c 
csc w + w' N + N'' 



(26) 



where A c is the area of the cavity. It is assumed that 
i osc > Lb, where 



ttA c /P c 



(27) 
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FIG. 1: Pair of orbits with a loop. (Neglecting the crossing 
region, we distinguish between the tails tl, t2 and the loop 



Spin-dependent quantum corrections to 
transmission and reflection 



We now compute the spin modulation factor for the 
leading contributions to the energy-smoothed T and 1Z, 
identified in Sec. lIIIBl First, the case with time-reversal 
symmetry<2& is considered: 

(i) For the classical part we find, using the unitarity 
of K 7 , that the modulation factor A4 7i7 = Tr (K^Kj!) = 
2s + 1 reduces to the trivial spin degeneracy. 

(ii) For the diagonal quantum correction the result is 



is the average distance between two consecutive bounces 
at the boundaries^ and P c is the perimeter. In an ar- 
bitrary billiard the last expression is true if the average 
is taken over the ensemble of chords with random ini- 
tial position and boundary component of the velocity. In 
ergodic billiards the average can, alternatively, be calcu- 
lated along almost any trajectory. 

(ii) The diagonal quantum correction is defined for re- 
flection only. It contains the terms with n = m and 
7' = 7 _1 , where 7 -1 is the time reversal of 7 (Ref. Esh ^i 
Clearly, when n and m are different channels, the or- 
bits 7(71,771) and 7'(n, m) cannot be the mutual time- 
reversals, since reversing the time would exchange n 
and rh. Again, the two actions are equal, and the re- 
sult for an ergodic system without spin reads^ 



sn 



(0) 



N 



ag N + N 1 



(28) 



(iii) The loop contribution consists of pairs of long orbits 
that stay close to each other in the configuration space. 
One orbit of the pair has a self-crossing with a small 
crossing angle e, thus forming a loop, its counterpart has 
an anticrossing, 85 Away from the crossing region the or- 
bits are located exponentially close to each other: they 
are related by time reversal along the loop and coincide 
along the tails^Sl (Fig. 0. The action difference for 
these orbits is of second order in e. For spinless chaotic 
systems with hyperbolic dynamics the loop terms in i|21|) 
yield^ 



NN' 



{N + N' 



sn 



(0) 
loop 



N{N+1) 
(N + N') 2 



(29) 



From here on, we will work in the limit N, N' 1. 
In this semiclassical regime (in the leads) the classical 
contribution 1)23(1 (of the order N) is much greater than 
the quantum corrections l|28|) and l|29l) (of the order N°), 
while higher-order loop corrections (of order ./V -1 and 
smaller) can be neglected. We note that the normaliza- 
tion is preserved order by order: 



(0) 
loop 



ST. 



(0) 



loop 



0(N~ 



(30) 
(31) 



Mn 



Tr(i^) 



(32) 



Kj, which follows 



It was taken into account that K' 

i' 

from the relation C 7 -i(t) = — C 7 (T 7 — t) and Eq. IJASJ. 

(iii) Assuming that the trajectories forming a loop pair 
(Fig. [JJ coincide along the tails £1, t2 and are mutu- 
ally time-reversed along the loop I, thereby neglecting 
the crossing region, we can represent the propagators as 
Kj = K t2 KiKti and -Ky = Kt2Ki~ 1 Kti- Hence, the 
modulation factor 



M 7 , y =Tv(Kf). 



(33) 



is independent of the tails. 



In the presence of a magnetic field the time-reversal 
symmetry is broken, and the preceding results ought to 
be adjusted. In this paper we consider a constant, uni- 
form, arbitrarily directed magnetic field B. Its compo- 
nent B z normal to the cavity is assumed to be weak 
enough, 87 so as not to change the classical trajectories 
in l)2ip. but only modify the action difference by the 
Aharonov-Bohm (AB) phase. We define the AB mod- 
ulation factor 



92 7 ,y = exp 



: A(5 7 -Sy) 



,-xi, 1 t-^-i ) , (34) 



Here, for a pair of trajectories 7 and 7' from the diagonal 
(loop) contribution, A = f A • d\/B z is the effective 
enclosed area, where the integral of the vector poten- 
tial A is taken along 7 (its loop part I), and <&o = hc/e 
is the flux quantum. 

With the Zeeman interaction included, Eqs. H32|) 
and l|33|) for the spin modulation factor are no longer 
valid. In fact, if we distinguish the SO and Zeeman terms 
in the effective magnetic field C 7 (t) = C^°(t) + C%(t), 
the diagonal (and, correspondingly, the loop) modulation 
factor can be written in the form 



M 



7i7" 



l(B) = Tr(fsT 7 if 7 ) 



(35) 



where 7 is a fictitious trajectory producing the field 
C 7 (i) = C^°(t) - C£(i). Clearly, in the absence of SO 
coupling, we have M^^-i(B) = 2s + 1, i.e., the Zeeman 
field alone does not affect the modulation factor. 
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In Refs. |48| and|49j the quantum corrections to trans- 
mission and reflection in the presence of magnetic field 
were calculated for an ergodic system. We extend their 
approach to a system with SO interaction. To this end, 
we consider the generalized modulation factor 



(36) 



The diagonal contribution can be computed using the 
sum rule ((24(1 . First, one averages (A4(p) ln > for the 
time-reversed pairs of trajectories and loops of a given 
length L. Thus, the average is performed over the en- 
semble of almost closed orbits. This restriction proves 
very important, since the average modulation factor for 
closed and open trajectories is different (see Sec|VJ. The 
average modulation factor M.Lp (L; B) is then further 
weighted with the length distribution ((25(1 . It can be 
shown that for the loop contribution in a hyperbolic 
system with a single Lyapunov exponent holds effec- 
tively the same procedure (see Ref. W\ and Appendix lB)l . 
Hence, the diagonal and the loop relative quantum cor- 
rections to reflection and transmission are equal and 
given by: 

JTZdtag (B) _ 6Tl loop (B) _ <57I oop (B) 



6K 



(o) 



5K 



(0) 
loop 



5%, 



(0) 



loop 



= (M<p(p)) L 

= - — / dL e - L / L ^ Mif (L; B). 

-^esc JO 



(37) 



The normalization condition l|18|l is preserved due to 
Eq. 

When the SO and Zeeman interactions are absent, the 
average modulation factor (2s + 1)^(L;B) in a chaotic 
system can be analytically estimated using the Gaussian 
distribution of enclosed areas^ 



P A (A;L) 



v/27rA2 L/Li 



exp 



2Ag L/Li 



(38) 



It depends on a system-specific parameter A , a typical 
area enclosed by an orbit during one circulation. This 
distribution does not depend on the incoming and out- 
going channel numbers and is valid for both closed and 
open trajectories. The average (|37[) yields^S*^ 



(M<p(B))l = 



2.s+ 1 



1 + B 2 L CSC /L b 



(39) 



where B = 2\/2itB z Ao/&o- Thus the quantum correc- 
tions have a Lorentzian dependence on the magnetic field. 
The result (|39|l is specific to chaotic systems as it depends 
on L csc — such a parameter is not relevant to extended 
disordered systems, while for regular billiards Pl{L) (in- 
troduced earlier) is usually a power lawa&. The increase 
of reflection (decrease of transmission) for B z = con- 
stitutes the effect of weak localization. A magnetic field 
destroys the time-reversal symmetry and, thereby, the in- 
terference between the mutually time-reversed and loop 
pairs of paths, thus diminishing the quantum corrections. 



The SO interaction may turn the constructive interfer- 
ence between the orbit pairs into destructive one. Since 
the sign of the quantum corrections in this case would 
be reversed, one speaks of weak antilocalization. In the 
following sections we study the transition from WL to 
WA and the related question of spin relaxation. 



IV. SPIN RELAXATION 

A. General discussion 

Equation 1)37(1 demonstrates that the modulation fac- 
tor A4tp (L; B) is a key to calculating the quantum cor- 
rections to the conductance. Hence, we will first examine 
M(p in detail. As a function of length, this quantity con- 
tains information about the average spin evolution along 
the trajectories of the system. Here, by the spin evolu- 
tion along a trajectory 7 we mean the change of the spin 
propagator K 7 (t). According to Appendix 1X1 it can be 
written in the form [Eq. IA6I1 ] 



K^t) = exp[-is • 77 (t)], 



(40) 



and, thus, depends on three real parameters: the rota- 
tion angle and the rotation axis given by the unit 
vector m 7 (t) = r\ (t)/r]j(t). Alternatively, one can pa- 
rameterize K~i using the elements of the corresponding 
SU(2) matrix [Eq. (U] 



(41) 



W (t) - p- iCT -'M t )/ 2 - ( 

w iW- e - y -b*(t) a ;{t) 

which are restricted by the condition det W 1 = |a 7 | 2 + 
|6 7 | 2 = 1. The two parameterizations are related by 
Eq. (|A4() . Clearly, W-y is the matrix representation of 
for spin s = 1/2. Instead of W 7 (t), we can consider the 

evolution of the spinor ip^(t) = (a 7 (t), — b* (t)j , start- 
ing from the spin-up state ip 7 (0) — (1,0) T . It is char- 
acterized by the spin direction n 7 (t) = [tp 7 (t)] T (T ipyif) 
and the overall phase (tr is the vector of Pauli matri- 
ces). For spin s > 1/2 these are the direction and 
the phase of a coherent state (see Sec. [nj. Note that 
n 7 (i) results from the rotation of n 7 (0) = (0,0,1) by 
the angle rjj(t) about m 7 (i). Sometimes it is conve- 
nient to represent W-y(t) by a trajectory on the three- 
dimensional unit sphere S 3 . For this purpose we define 
a four-dimensional unit vector 



£ 7 (i) = (— Im& 7 (i), — Rc6 7 (<), — Im a 7 (t), Re a 7 (i)) 



m 7 (i) sin 



(42) 



The trajectory starts at the "north pole," i.e., £ 7 (0) = 
(0,0,0,1). 

Using the propagator matrix 141(1 , the modulation fac- 
tor Eqs. ((32(1 and l(33|) for spin 1/2 can be expressed as 



M = Tr (W 2 ) = 4£ 4 2 - 2 = 2cost? 



(43) 
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where £4 is the fourth component of £ (in Appendix O an 
arbitrary s is considered). To simplify the notation, we 
dropped the subscripts labeling the trajectory, and the 
time t is the trajectory or loop time. 

For long orbits one expects that the spin state be- 
comes completely randomized due to SO interaction, if 
the particle motion is irregular. This means that all 
points £ € S 3 are equally probable, and, on average, 
£ 4 2 = 1/4 in the limit L — > 00. Hence, for B = the 
modulation factor A4(L) = A4(p(L;0) changes with L 
from the positive value A4(0) — 2 to the negative asymp- 
totic value M(oo) — —1 (cf. Ref . l68l) . The stronger the 
SO interaction, the shorter the length scale Lj^t of this 
change. If L esc <C L^, i.e., the particle quickly leaves 
the cavity, or the SO interaction is weak, then the rel- 
ative quantum corrections (|37|l are positive, giving rise 
to WL. In the opposite case of strong SO interaction or 
long dwell times, the relative quantum corrections are 
negative, leading to WA. For an arbitrary spin, Ai(L) 
changes from M(0) = 2s + 1 to M (00) = (-l) 2s (Ap- 
pendix Q. Thus, WA cannot be observed for an integer 
spin, at least, for L csc ^> Lm \M.{L) can, in principle, 
become negative at intermediate lengths]. 

For the rest of the paper we will consider the physically 
most important case of spin s = 1/2. If the Zeeman 
interaction is included, then Eq. (|35|1 yields 

M(B) =4&&-2£-£, (44) 

where £ belongs to the fictitious trajectory 7. In the ab- 
sence of Zeeman coupling, the vectors £ and £ coincide. 
Then the negative second term in Eq. I|44|) is responsi- 
ble for the WA, if the first term is, on average, small 
enough due to SO interaction. An admixture of a mod- 
erate Zeeman coupling destroys the correlation between £ 

and thereby reducing the average product £ • Thus, 
an external magnetic field suppresses WA in two ways: 
the AB flux breaks down the constructive interference 
between the orbital phases and the Zeeman interaction 
affects the spin modulation factor. As we know, the for- 
mer mechanism inhibits the WL, as well. 

The spin propagator K^(t) can be used not only in the 
calculation of the quantum corrections (|37|l — it also pro- 
vides information about the spin relaxation along clas- 
sical trajectories, which is of separate interest. The re- 
laxation of the spin direction can be described by the 
z component of the vector n 7 : 

n 2 = 2(£ 2 +£ 4 2 )-l. (45) 

The ensemble average 777 (L) varies from 777 (0) = 1 to 
777" (00) = 0, if the memory of the initial spin direction is 
completely lost for long orbits. The typical length scale 
of this decay can be different from Lm , because M. (L) 
depends on the phase of the spin state, as well as on 
its direction. Moreover, the length scale of 777 (L), as 
defined by Eq. I|45|l , depends on the choice of the quanti- 
zation axis. An invariant measure of the spin relaxation 



is given by £4 (L) or, equivalently, by A4(L). The differ- 
ent relaxation rates of 777/ (L) and A4(L) are observed in 
two-dimensional systems with the Rashba and the Dres- 
selhaus SO coupling (see Sees. II V 51 and Ivjl . 

B. Example: diffusive systems 

In three-dimensional extended diffusive conductors 88 
the directions of the effective magnetic field C 7 (i) be- 
fore and after a scattering event can be assumed uncor- 
rected 2 ^. We model this by keeping |C| = const and 
changing the direction of C randomly at identical time 
intervals equal to the elastic scattering time r. The spin 
propagator for the jth time interval is Kj = exp (—is • 
nij |C| r), j = 1,2,..., where mj is a random unit vec- 
tor. The position on S 3 after the first time interval is, 
according to Eq. (|4*2")l . £(t) = (m.i sin— ij-^, cos — • 

Thus, a trajectory on the sphere, starting at the "north 
pole," traverses an arc of length |C| r/2 along a randomly 
chosen great circle. During the second time interval, the 
trajectory starts at £(r) and moves along another random 
great-circle segment, and so on. Clearly, £(t) follows a 
random walk on S 3 . In the continuous limit |C| r <C 2-7T, 
its probability density satisfies a diffusion equation. Solv- 
ing this equation (Appendix [D| we find that the average 
modulation factor for trajectories of time i, 

M diff , 3D (t) = 3e-^ c l 2rt -l, (46) 
and the average spin polarization, 

Md iff ,3DW = ^ |C|2Tt , (47) 

exhibit the same relaxation rate. Note that Eq. (|37fl is 
not valid in diffusive systems. The modulation factor l|46() 
is equivalent to the result of Eq. (10.12) of Ref. l29i 

In two-dimensional diffusive systems with Rashba or 
Dresselhaus interaction it is reasonable to assume that C 
acquires a random direction in a two-dimensional plane. 
In this case the walk on S 3 is not fully random. As our 
numerical simulations show (Fig.[2J), the modulation fac- 
tor is reasonably well described by Eq. H46|) . However, 
the off-plane polarization 777 (t) relaxes faster in two di- 
mensions [cf. Eqs. (34) and (35) of Ref. This is not 
surprising, since, obviously, £3(1") = in 2D, but not 
in 3D. 



V. RASHBA AND DRESSELHAUS 
INTERACTION: SPIN RELAXATION 

A. Effective magnetic field 

We apply the general theory of the previous sections 
to ballistic quantum dots with Rashba 3 and Dresselhaus 4 
SO interaction. Both contributions are usually present 
in GaAs/AlGaAs heterostructures. Their strength ratio 
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FIG. 2: Average spin modulation factor A4(t) and spin pro- 
jection nZ (t) in diffusive systems. The effective magnetic 
field C(t) changes its direction randomly at equal time inter- 
vals r. Its magnitude is kept constant and is equal to 0.3 /r in 
this example. The analytical expressions 14611 and 147H (solid 
curves) are compared with the results of numerical simula- 
tions for C in two (dash-dotted curves) and three (circles) 
dimensions. The average was performed over 10 5 random se- 
quences ("trajectories"). 

can be experimentally varied, e.g., by tuning the Rashba 
SO strength through an additional gate voltage 70 , When 
the two-dimensional electron gas lies in the (001) plane 
of a zinc-blend lattice, the effective magnetic field C = 
Cr + Co in the Hamiltonian consists of 

Cfl = ^vxe z , (48) 

A.R 

C D = - — (v x e x - v v e y ), (49) 
Ad 

where the x and y axes are chosen along the [100] 
and [010] crystallographic directions, respectively, and 
v — (p — eA/c)/M is the (Fermi-) velocity operator de- 
pending on the effective mass M. In Eq. [Eq. (|4Tj|l ] 
the Rashba (Dresselhaus) interaction, usually character- 
ized by the constant an (an), is measured in terms of 
the inverse spin-precession length A^^-j = ctR^M/irK 2 . 
In billiards the natural dimensionless parameter is 

6 R[D) =2itL b /k R{D) . (50) 

It signifies the mean spin-precession angle per bounce if 
only one type of SO interaction is present. 

As can be seen from Eq. I|48|l . the effective Rashba 
magnetic field C R (t) generated by a particular trajectory 
points perpendicular to the velocity v(t). The directions 
of the Dresselhaus field Cn(t) and v(i) are symmetric 
with respect to the x axis [Eq. ljl9"|l ]. Hence, C R (t) and 
Co(t) always point symmetrically with respect to the 
[110] direction, labeled here by X (Fig. [3J. As a con- 
sequence, the total field C(t) is reflected about X un- 
der the exchange A R ^ A D in Eqs. and (EE9"|) . This 
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R ^X[U0] 

FIG. 3: The Rashba effective field Cr is normal to the veloc- 
ity v, while the directions of the Dresselhaus field Cd and v 
are symmetric with respect to [100]. Thus, the directions of 
Cfl and Cd are symmetric with respect to [110]. 

means that the modulation factor M. (t) and the polariza- 
tion projections n z (t) and nx(t) are preserved under this 
transformation. For example, systems with only Rashba 
or only Dresselhaus interaction have identical spin evo- 
lution, if the coupling strengths are the same. 

It is sometimes convenient to work in the coordinate 
frame of X and Y = [110] (cf. Ref. Et). The projections 
of the effective magnetic field on these axes are given by 

C x = 2wv Y /A x , C Y = 2ttv x /A y , (51) 

where 

A^^A^ + A^ 1 , A^A^-A^ 1 (52) 

are the effective inverse precession lengths. As above, 

we can define dimensionless parameters 9 x /y) = 
2irL b / A X (y)- 

B. Numerical study 

The computation of the spin evolution in billiard cav- 
ities is relatively straightforward, since the classical tra- 
jectories there are sequences of straight segments. If 
only the uniform Rashba and Dresselhaus interaction is 
present, the effective magnetic field Cj is constant along 
the jth segment (the segment velocity Vj is constant, 
moreover, its magnitude v is the same for all j due to the 
energy conservation). The spin-propagator matrix 141|) 
for a trajectory, 

w = e -io- v /2 = Wl ..-Wi, (53) 
is a product of the respective matrices 

Wj = e- iCT -V 2 = e -i°"<W2 = 1, . . . , (54) 

for the I orbit segments. In practice, it is convenient to 
remove the velocity dependence in Eqs. (|51|l by using the 
displacement Ar\, = AX,,- e x + AY,- ey = Vj tj instead of 
the segment time tj . Thereby the rotation vector can be 
expressed as 

/AY,- AX, \ . , 
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FIG. 4: Billiard geometries: desymmetrized Sinai (DS) bil- 
liard, desymmetrized diamond (DD) billiard, desymmetrized 
Bunimovich (DB) stadium billiard, quarter circle (QC), rect- 
angle, and circle. The leads are numbered for future reference. 



It follows from this equation that rescaling of the system 
size and the spin-precession lengths by the same factor 
does not change the spin relaxation. In other words, 
given the shape of the billiard, the averages and M. 
as functions of L /Lb (computed below) depend only on 
the angles Or and 0r>- 

We performed a systematic numerical study of spin 
relaxation for several billiard geometries (Fig. 0} repre- 
sentative for systems with chaotic and integrable classi- 
cal dynamics. The desymmetrized Sinai (DS) billiard, 
the desymmetrized diamond— (DD) billiard, and the 
desymmetrized Bunimovich (DB) stadium billiard repre- 
sent chaotic cavities. The quarter circle (QC), rectangle, 
and circle are integrable. The average spin relaxation is 
computed for the closed versions of these billiards. 

Figures |S] and depict the average spin relaxation de- 
scribed by n^(L) and its invariant counterpart, JA(L), 
for the chaotic DS and integrable QC billiard. The av- 
erage is performed over ensembles of open (Fig. [5J and 
closed (Fig. |BJ) trajectories of length L starting at ran- 
dom position at the boundary with a random boundary 
component of the velocity. The strength of the Rashba 
interaction is chosen as 0r/2tt — 0.2, and the Dressel- 
haus interaction is absent (or, equivalently, 9d/2t: = 0.2 
and Or = 0). 

In Fig. [SI the numerical results for an extended two- 
dimensional diffusive system, where Lb is identified with 
the elastic mean free path vt, are shown for comparison. 
We observe that on the scale of L ~ Lb the spin relax- 
ation is the same in all three examples. Indeed, before the 
first collision with the boundary or a scatterer, the parti- 
cle moves along a straight line, irrespective of the system 
it belongs to. On longer length scales relaxation in an ex- 
tended diffusive system is much stronger than in confined 
systems. Moreover, in the integrable billiard saturation 
takes placed We also note that n^(L) in the chaotic bil- 
liard, similarly to the diffusive system (Sec. lIVB|) . relaxes 
to its asymptotic value faster than A4(L). 




FIG. 5: Average spin projection n7(L) and modulation fac- 
tor A4(L) for the closed chaotic DS (solid curves) and inte- 
grable QC (dashed curves) billiard (see Fig. |3J . Each data 
point represents the average over 50,000 open trajectories 
with random initial values at the boundary. The numerical 
results for a two-dimensional extended diffusive system (dash- 
dotted curves), where Lf, is identified with the elastic mean 
free path vt, are shown for comparison. The SO-coupling 
strength is 6r/2tt = 0.2 and 6n = 0. 
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b 

FIG. 6: Average spin projection nJ(L) (dotted curves) and 
rescaled modulation factor [A4(L) + 1]/3 (solid curves) for the 
closed chaotic DS and integrable QC billiard. Each data point 
represents the average over 5,000 closed trajectories started 
at random at the boundary. The SO-coupling strength is 
9 R /2iv = 0.2 and 8 D = 0. 
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FIG. 7: The average modulation factor M(V) for the closed 
chaotic DS, DD, and DB billiards. In the latter case the 
ratio of the upper straight side to the radius was taken to 
be 1 (DB, 1) and 3 (DB, 3). Each data point represents 
the average over 50,000 open trajectories started at random 
at the boundary. The SO-coupling strength is 8r/2tv = 0.1 
(four upper curves) and 0.2 (four lower curves), while 6n = 0. 



The ensemble of closed orbits (Fig. is responsible 
for the quantum corrections to transmission and reflec- 
tion (Sec. IIIip . We find that the relaxation in this case is 
much slower than for the ensemble of arbitrary trajecto- 
ries. Remarkably, the spin projection nj and the rescaled 
modulation factor (~M + l)/3 are hardly distinguishable. 

The spin evolution in several chaotic systems is com- 
pared in Figs. [7| and 03 for the ensembles of open and 
closed trajectories, respectively. All the billiards show a 
qualitatively similar behavior. The DS and DD billiards 
have about the same relaxation rate. In the DB bil- 
liards the relaxation rate grows continuously, starting 
from zero, as the ratio of the upper straight side to the 
radius increases. 

In Fig.OHthe modulation factor averaged over the open 
trajectories is presented for the integrable QC and rect- 
angle billiard. Both systems are characterized by the 
saturation of spin relaxation and persistent long-time os- 
cillations. The saturation level decreases down to —1 as 
the SO coupling becomes stronger. The circle billiard, 
which shows a non-typical relaxation pattern, is consid- 
ered in Sec. WU 

When the Rashba and the Dresselhaus couplings work 
simultaneously, they mutually counteract their effects on 
the spin relaxation 90 (Fig. I10|l . In the extreme case Ar = 
Ad, i.e., Ay 1 = 0, the effective field C is always parallel 
to the X axis2i (Fig.[3|). Hence, the propagator matrices 
in Eq. (|53|l commute, and the rotation vector becomes 

/ 

77 = Vj = 27t {AY/ A x )e x , (56) 
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FIG. 8: Same as in Fig. but for the ensemble of closed 
orbits started at random at the boundary. Each data point 
represents the average over 5,000 trajectories. 
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FIG. 9: The average modulation factor M(L) for the closed 
integrable QC (solid curves) and rectangle (dashed curves) 
billiards. Each data point represents the average over 50,000 
open trajectories started at random at the boundary. The 
SO-coupling strength is 9r/2tt — 0.1 and 0.3, while 6d = 0. 

where AY — Ylj=i is the Y displacement for the 
trajectory. According to Eq. H43|) . the modulation factor 
is then 

M = 2 cos (2tt \AY\/A x ). (57) 

On the long- length scale L Lb, \AY\ varies from orbit 
to orbit between and the system size. Clearly, the aver- 
age M. should be independent of the orbit length L. This 
explains the saturation in Fig. ^| The saturation level 
decreases from 2 to as 9r changes from to oo. For 
closed orbits we have AY" = 0, and, therefore, Ai = 2. 
If the spin is initially polarized in the X direction, the 
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FIG. 10: The average modulation factor A4(L) for the closed 
DS billiard at different strengths of the Rashba and Dressel- 
haus interaction. Each data point represents the average over 
50,000 open (dashed curves) or 5,000 closed (solid curves) tra- 
jectories started at random at the boundary. The values of 
(0r/2tt, 9d/'2h) are shown on the graph. The leads are paral- 
lel to the [100] direction. Note that M — 2 for closed orbits, 
when Or — On- 



polarization does not change with L, i.e., nx — 1. On 
the other hand, Eqs. JUJl and (gSJ yield n z = M/2 if the 
trajectory starts with n z = 1. This demonstrates that 
the spin-relaxation measure n^(L) depends on the choice 
of the quantization axis, as was mentioned in Sec. II VI 



C. Limit of slow spin motion 

Above we presented a number of numerical obser- 
vations regarding the average spin relaxation in two- 
dimensional billiards. A further insight into the connec- 
tion of the spin dynamics to the characteristics of the or- 
bital motion can be gained within the limit of slow spin. 
In this approximation the period of spin precession is, by 
definition, much longer than the time scale on which the 
orbital momentum changes. In billiards this requirement 
takes the form \C\L^/v <C 27r, or, 9xy *C 2-7T. 



1. Rotation-angle expansion 

The rotation angle rj [Eq. {53} ] contains all the in- 
formation about the spin evolution along a particular 
trajectory. It essentially depends on the rotation an- 
gles rjj for the straight segments of the trajectory. In 
turn, the angles rjj are directly related to the orbital dis- 
placements via Eq. I)55|l . Thus, a more explicit expression 
of 7] in terms of the rjj is desirable in order to establish 
the link between the geometry of orbital motion and the 
spin rotation. For this purpose we employ the Baker- 
Campbell-Hausdorff (BCH) formulaSS for the product 



of the exponentials of two matrices (or operators), 

exp (P) exp (Q) = exp ^ rA , (58) 

where Ri's are homogeneous polynomials of degree i in 
P and Q. The first three of them are given by 



Ri 
R 2 

R 3 



P 
1 

2 
1 

12 



Q, 
[P,Q], 

[P-Q,[P,Q]] 



(59) 
(60) 

(61) 



([,] denotes a commutator). The BCH formula can be 
used to calculate the product of the segment propaga- 
tors Wj in Eq. i|53|) . However, only in the limit of slow 
spin, when \r\j\ -C 1, the first few contributions Ri are 
sufficient. The rotation angle for a trajectory, 



T) = »7 (0) + St] 1 - + 5T7 11 , 
comprises three parts. The lowest-order term, 



V 



(0) 



= 9tt 



V A 



ex 



x 



AX 



e Y 



(62) 



(63) 



is a vector sum of the segment rotation angles (AX = 
5Z 7 -=i &Xj). The correction normal to the billiard plane, 



dr) 1 - 



A X A Y 



Ae z 



0(Ax 4 Y ), 



(64) 



is proportional to the effective enclosed area A (Ap- 
pendix El ■ For closed orbits, A coincides with the area 
defined below Eq. I|34|) for a returning orbit or a loop. In 
an open trajectory, A is the area of its "closure" obtained 
by connecting the endpoints with a straight line. The in- 
plane correction Srj^ is of the order A^ Y - I n general, the 
in-plane (normal) contribution to r] contains odd (even) 
powers of r\j . 

Expressions (|63|l and (|64|l help to interpret some of the 
results of Sec. EH 

(i) The suppression of spin relaxation along open tra- 
jectories in confined systems, as compared to extended 
diffusive systems (Fig. |SJ|, can be deduced from the ex- 
pansions 



M/2 



1 



(0) 



0(Axy), 



(65) 



that follow from Eqs. (021 an d In billiards |tj(°) | in- 
creases from zero to its saturation value of the order 
of {Ox.y) 2 on the length scale of the system size. The 



X.Y 



further relaxation on the scale L > Lb is due to the A 
order terms. In diffusive systems, on the other hand, 

r/(°) grows linearly with length. 
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(ii) For the closed orbits one sets r/°) = 0. Hence, the 
short-scale relaxation does not show up in the ensemble 
averages M{L) and nj(L) (Fig.|fjJ). If the in-plane contri- 
bution Srj" is neglected, i.e., \Sr]"\ <C 1, the modulation 
factor becomes 



M ~ 2cos|<5r7- L |, 
and the generalized modulation factor is 



Mif ~ exp 



V $0 A x A y 



(66) 



(67) 



In a chaotic system, averaging with the area distribu- 
tion (SU yields 



Mtp(L;B z ) ~ J^exp 



-(S±A- 2 ) 2 ^ 



(68) 



where A -2 = 2\/27r 2 Ao/AxAy. Clearly, the normal con- 
tribution Sr) 1 - alone cannot make Aiip negative. For suf- 
ficiently large lengths L, the components Sr) 1 - and Srj^ 
become comparable and reverse the sign of the modu- 
lation factor. With the help of Eqs. 03) and (@2J), we 
obtain the spin polarization 



n z = 1 — 2 (m 2 + m 2 



1 — — — I Sin 



|<5r, 



-LI 2 



(69) 



The Sry" component is responsible for the rotation of 
an initially spin-up state. This means that for suffi- 
ciently weak SO interaction and short L, relaxes 
slower than M. (Fig. 111(1 : the relaxation rates are of the 
order of A_^y and A^ Y , respectively. For stronger in- 
teraction the difference is not noticeable (Fig. EJ) ■ 

(iii) In integrable billiards the strong area cancellation 48 
along trajectories contributes to the saturation of spin 
relaxation. The circle billiard makes an exception: here 
all the trajectories accumulate area linearly in time. For 
an orbit having the shortest distance r from the center, 
the enclosed area is, on average, A ~ ±rL/2 (the signs 
denote the two senses of rotation). Averaging the mod- 
ulation factor H67f) over all ±r, i.e., assuming that all 
angular momenta are equally possible, we find 4 ^ 



sma± 
a± 



a± = 2irRL ( — ^ ± 



$o A x A y 



(70) 



where R is the radius of the circle. [Note that a sim- 
ilar average yields Lb = wR/2, in accordance with 
Eq. J23-] Surprisingly, for a sufficiently small SO cou- 
pling, the average over open trajectories agrees better 
with Eq. H70[). than the average over the closed orbits 




10000 



FIG. 11: Average spin projection nj(-L) (dashed curves) and 
rescaled modulation factor [j\4(L) + l]/3 (solid curves) for 
the closed chaotic DS billiard. Each data point represents the 
average over 15,000 closed trajectories started at random at 
the boundary. The SO-coupling strength 9r/2tt is shown in 
the graph and 8d = 0. 



started at the boundary [Fig.^](a)]. This happens be- 
cause the J] contribution is small, but different angu- 
lar momenta are not equally represented in the ensem- 
ble of closed orbits. As well as Eq. l|55)l. Eq. {7UJ| docs 
not describe the full relaxation at large SO interaction 

[Fig. in (b)]. 



2. Unitary transformation of the Hamiltonian 

It is possible to transform the Rashba-Dresselhaus 
Hamiltonian [Eqs. © and l|51|l ] to the form in which 
the SO interaction is weaker by a factor of &x,y -C 1 
fRefs. Ell and 113) . Considering the case B = 0, we start 
with the Hamiltonian 



H 



2M 



A 



x 



A 



Y 



and apply the unitary transformation 



V = exp 



A 



x 



A, 



(71) 



(72) 



to it. In the limit of slow spin motion, the exponent 
can be expanded in powers of A^y. Keeping only the 
terms quadratic in 9x,Y and linear in H/Ax,y [due to 
the weak-coupling assumption i|l(J|) ]. we obtain the new 
Hamiltonian 



H = V^HV 



P 

2M 



- (aw 



A (2) 



o 



'X,Y, V A 



M 

h 



X.Y , 



(73) 
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FIG. 12: The average modulation factor A4(L) for the closed 
circle billiard of radius 7?. Each data point represents the 
average over 50,000 open (solid curves) or 5,000 closed (dot- 
ted curves) trajectories started at random at the bound- 
ary. The results are compared with the analytical expres- 
sion I70H (dashed curves). The SO-coupling strength is 
e R /2n = ttR/2A r = 0.1 (a) and 0.3 (b), while 9 D = 0. 



It contains the spin-dependent "vector potentials" 

A« = -A ( e zXr)ffl , (74) 

Ax Ay 

A (2)_ 4 ^ (c xr) (Z2L_XZL\ (75) 
A - 3A x A Y [ezXr) { Ay A X ) (?5) 

of the order of (h/Ax.y) &x,y and (h/Ax,y) #x y , respec- 
tively. 

As any two Hamiltonians connected by a unitary trans- 
formation, H and H must yield the same total transmis- 
sion T and reflection TZ. While H and H have identical 
orbital parts, the SO interaction in the latter is much 
weaker. It is important, therefore, to understand the 
role of the transformation V within the framework of the 
semiclassical Landauer formula developed in the previous 
sections. Working in the weak-coupling regime (|10l) and 
treating the orbital motion semiclassically, we consider 
the unitary transformation 



V(t) = exp < — in 



Y(t)a x , X{t)a Y 



A 



x 



A^ 



(76) 



acting in the spin space. It is generated by a classical 
trajectory [r(f),v(f)]. The spin-propagator matrix 141|) 
transforms as 



[This property follows from the spinor transformation 
$(t) = V^{t)ip{t) and the definition ip(t) = W(t)$(0).] 
For a closed orbit of time t, we have V(t) = V(0). Hence, 
the modulation factor for a pair of trajectories from the 
diagonal or the loop contribution (Sec. IIIIL 



M(t) = Tr W(t) =M(t) 



(78) 



is preserved under the transformation. Thus, H and 

H return the same conductance also in the semiclassi- 
cal approximation^ This result demonstrates the sig- 
nificance of closed trajectories for the transmission and 
reflection. 

According to Eq. (|78|l . the modulation factor for closed 
orbits can be determined directly from the semiclassical 

version of H . Neglecting A^ 2 ^ , we compute the propaga- 
tor as 



W(t) ~ exp 



exp 



- / AM-vdf 
2tt 2 . 



Ax Ay 



(79) 



Then Eq. ijrjU) for the rotation angle and Eq. i|(77|) for the 
modulation factor (with B = 0) follow. 

For the open trajectories, the leading-order (in A^ Y ) 
contribution from Eq. (|77|l is 



W(t) ~ V{t)V*{0) 

~Y(t)-Y(0) 



exp 



Ax 



X(t) - X (0) 

Ay 



■cry 



(80) 



As a consequence, Eqs. (|63|) and (|65|l are obtained. 



W(t) = V\t)W(t)V(0). 



(77) 



VI. RASHBA AND DRESSELHAUS 
INTERACTION: QUANTUM CORRECTIONS TO 
TRANSMISSION AND REFLECTION 

A. Results of the semiclassical theory 

As was shown in the previous sections, the average 
relaxation of the spin modulation factor as a function of 
orbit length is closely related to the type of the classical 
dynamics of the system. In view of Eq. H37jl. the effect of 
the classical orbital motion is also felt in the magnitude 
and sign of the quantum correction to the conductance. 

In the examples below the following numerical pro- 
cedure is employed to calculate the relative quantum 
corrections: The trajectories are started at the bound- 
ary within a given lead. Their initial coordinate and 
velocity component at the lead cross section is cho- 
sen at random. If an orbit returns to the same lead, 
the generalized modulation factor l|36|) for this orbit 
and its time-reversed partner is recorded. The aver- 
age over these modulation factors yields the diagonal 
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FIG. 13: Relative quantum correction to reflection (or trans- 
mission) coefficients vs. Rashba-interaction strength 8a (8d = 
0). The data are shown for the chaotic (DS and DD) and inte- 
grable (QC, rectangle, and circle) billiards with P c /(w+w') — 
90. Each data point represents the average over 50,000 orbits. 
The details on the numerical procedure are given in the text. 
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FIG. 14: The relative quantum correction to reflection (or 
transmission) coefficients as a function of 8r for the chaotic 
DS billiard with P c / (w + w') — 90. Solid curves: fixed ratios 
0r/8d (shown in the graph). Dotted curve: fixed 9d = 7T. 
Each data point represents the average over 50,000 orbits. 



contribution^ <57?.diag (B) / SlZ diag , which, by virtue of 
Eq. I|37|) . is equal to the total relative quantum correc- 
tions Sn(B)/STZ {0) = 8T(B)/8T^ for a chaotic system. 
For non-chaotic billiards there exist no analytical ap- 
proach to the loop statistics. Hence, there we use the 
diagonal part as an estimate for STZ(B) / STZ^. 

In Fig. I13l the relative quantum correction is shown as 
a function of the Rashba-interaction strength 9r (while 
8d = 0) for several billiards. Here and in the following 
figures the results for lead "1" (Fig.^J are presented. All 
the billiards have the same ratio P c /(w + w'), which in 
chaotic systems gives the escape length in units of Lb- As 
9fi increases, the relative quantum corrections in chaotic 
cavities decrease and, eventually, change sign. The level 
STZ/SU^ = corresponds to the WL-WA transition. In 
the integrable QC and rectangle billiards no WA occurs, 
while the circle billiard shows a non-typical pattern. 

The interplay between Rashba and Dresselhaus inter- 
action in a chaotic billiard is investigated in Fig. El As 
was already noted in Sec. IV Bl the two interactions tend 
to compensate each other (thus suppressing the WA) as 
their ratio gets closer to unity. When Or = 9 m there 
is no manifestation of SO interaction in the conductance 
[57?.(B = 0)/STZ^ ~ 2], since, effectively, no spin relax- 
ation is associated with closed trajectories. 

Owing to the scaling property mentioned below 
Eq. (JS5J, the quantum corrections for a given billiard 
shape depend only on 9r and 9b- These parameters are 
proportional to the size of the system. Thus, in a material 
with definite physical SO-coupling constants cxr and am 
Fig- E] (curves with constant 9r/9d) predicts the sup- 
pression of WA with the decreasing size of a chaotic quan- 
tum dot (also observed in experiments^) . Preserving the 
shape of the cavity means, in particular, that P c /(io+u/) 



is kept constant when the dot size is changed. Alterna- 
tively, one can fix w + w': this leads to a disproportionate 
decrease of L esc when the size is reduced, thereby sup- 
pressing WA even stronger. 

In experiments, for fixed 9 m 9r can be tuned through 
an additional back-gate potential. The dotted curve in 
Fig. El shows the result of such a parameter variation, 
including the associated WA- WL-WA crossover. 

The effect of Zeeman interaction in the chaotic DS bil- 
liard is considered in Fig. 1151 The interaction strength is 
measured by the spin-precession angle per bounce 9z = 
Lfc|Cz|/|v|, where Cz is the Zeeman contribution to the 
effective field C. We observe a complete suppression of 
WA even by weak Zeeman coupling. In general, the 
9z dependence of SIZ/SIZ^ is very complex. It shows 
strong anisotropy with respect to the direction of the 
applied in-plane magnetic field, especially, for weak SO 
interaction. Note that the symmetry of the spin relax- 
ation under the Ar <-> Ko exchange (Sec. IV Af) is lifted, 
unless Cz lies in the (110) plane. 



B. Comparison with full quantum mechanics 

The semiclassical results are compared with quan- 
tum calculations in Figs. El an d El The complete 
quantum transmission and reflection amplitudes t na i, ma 
and r na i_ mcr are computed numerically by using a re- 
cursive Green-function technique (see Ref. and ref- 
erences therein for details). This approach is based on 
a tight-binding model arising from the real-space dis- 
cretization of the spin- dependent Schrodinger equation in 
a two-dimensional geometry. Here, the standard on-site 
and hopping energies present in a spinless tight-binding 
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FIG. 15: The relative quantum correction to reflection 
(or transmission) coefficients as a function of the Zeeman- 
coupling strength 8z for the chaotic DS billiard with P c /(w + 
w') — 90. The in-plane magnetic field is parallel (solid curves) 
or perpendicular (dotted curves) to the leads' direction (see 
Fig. 0J. Two Rashba-coupling strengths 6 R /2n = 0.1, 0.3 
(shown in the graph) are considered. The Dresselhaus inter- 
action is absent. Each data point represents the average over 
50,000 orbits. 
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FIG. 16: Relative quantum correction to reflection coefficient 
vs. Rashba-interaction strength 8r (8d = 0) in the DS billiard 
with P c / (w+w') — 30. The quantum reflection was computed 
for the energy corresponding to N — N' — 6 open channels in 
the leads. The relative quantum correction (solid curve) was 
obtained after estimating the classical reflection 1Z c i — 6.49 as 
described in the text. In the semiclassical case (dashed curve) 
each data point represents the average over 50,000 orbits. 



model are replaced by 2 x 2 matrices accounting for the 
spin degree of freedom^ In combination with a recur- 
sive algorithm, we obtain the Green function describing 
spin-dependent transport between the leads by solving 
the (implicit) Dyson equation. The amplitudes t na > tma 
and r na ) . ma are evaluated by projecting the Green func- 
tion onto a complete set of asymptotic states defining 
the incoming/outgoing spin and orbital channels in the 
leads. The method has the advantage of large flexibil- 
ity: the quantum amplitudes can be equally calculated 
for geometries of arbitrary shape. Moreover, disorder, 
magnetic fields and spin-dependent interactions can be 
easily introduced and modified. 

The dependence of STZ/STZ^ on the Rashba-coupling 
strength in the DS billiard is presented in Fig. ^] 
Compared to Fig. EH here we choose a smaller ratio 
Pc/(w + w ') m order to reduce the computing time in 
the quantum case. There is a good overall agreement 
between the two curves, especially, above the WL-WA 
transitional region. Since full quantum calculations give 
only the total reflection, one is faced with the problem 
of separating the quantum corrections from the classical 
part. Instead of evaluating 1Z C \ directly, we estimate it by 
requiring that the ratio 51Z{9r = 0)/5TZ(0r 3> 1) is the 
same in the semiclassical and quantum consideration. 

We used the estimated value of 1Z C \ for calculating 
the quantum corrections in the presence of magnetic 
flux B z A c /$o (Fig. I17|) . In this graph the quantum re- 
sults exhibit additional oscillations around the semiclas- 
sical curves, both with and without the SO interaction. 
The observed deviations from the semiclassics might be 
due to the fact that, for the energy range considered in 
the quantum calculations, the bending of classical trajec- 
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FIG. 17: The relative quantum correction to reflection coeffi- 
cient vs. magnetic flux in the DS billiard with P e /(tu + w') = 
30. The quantum reflection was computed for the energy cor- 
responding to N = N' = 6 open channels in the leads. In the 
semiclassical case each data point represents the average over 
50,000 orbits. The Rashba interaction was kept constant at 
0r/h — (solid curve: semiclassical / circles: quantum); 0.5 
(dash-dotted curve / triangles); 0.66 (dotted curve / squares); 
1.33 (dashed curve / diamonds). The Dresselhaus and Zee- 
man couplings were not included. 



tories caused by a Lorentz force limits the applicability of 
the theory in Scc. lIII 01 By the same token, the change of 
the trajectories resulting from the SO interaction might 
be responsible for the discrepancy in Fig. ^] Neverthe- 
less, the semiclassical techniques prove to be an effective 
tool for making reasonable predictions with a minimum 
of computational power^ 



16 



The maximum in 6TZ(B Z )/STZ {0 \ clearly noticeable at 
the intermediate SO interaction (Fig. I17fl . can be at- 
tributed to the cancellation between the external flux 
and the flux of the SO vector potential (|74|! 42 . Equiv- 
alently, the integration 1|37|) with the modulation fac- 
tor (S| yields [cf. Eq. J2jJ)] 



811(B) 



E 

± 



l + (B±A-^L csc /L b 



(81) 



Thus, the maxima are located at B~ ±A~ 2 . Note that 
the above expression, being positive, does not describe 
the WA, which appears due to the drj^ spin rotation (see 
Sec.lVCl 



the restriction to closed orbits is lifted, if the leading- 
order Rashba/Dresselhaus terms are removed from the 
Hamiltonian by a gauge transformation. 

The semiclassical results for a chaotic billiard with 
moderate spin-orbit interaction and weak magnetic field 
show a good agreement with the quantum calculations. 
To account for the deviations at larger values, it may be 
necessary to include the distortion of classical trajecto- 
ries in the semiclassical analysis. Possible extensions of 
the theory could be based on the extended-phase-space 
approach to spin-orbit coupling^iS. 
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VII. SUMMARY AND CONCLUSIONS 

The semiclassical Landauer formula with spin ex- 
presses quantum corrections to transmission and reflec- 
tion in ballistic quantum dots as a sum over pairs of 
classical trajectories or their loop parts, related by time 
reversal. The effect of the spin-orbit and Zecman interac- 
tion, as well as the Aharonov-Bohm phase, can be taken 
into account by averaging of a modulation factor over the 
ensemble of closed orbits. The change of classical trajec- 
tories by the spin-orbit interaction and the Lorentz force 
can be neglected, provided the energy of the particle is 
sufficiently high. With the increasing spin-orbit coupling 
strength the quantum corrections may reverse their signs, 
and weak localization will become weak antilocalization. 
Whether such a transition takes place, depends on the in- 
terplay between the average spin relaxation and the dwell 
time in the quantum dot. 

The spin relaxation, particularly for Rashba and Dres- 
selhaus interaction, is very sensitive to the character of 
classical dynamics, which is determined by the shape of 
the boundary. While in chaotic geometries spin eventu- 
ally completely relaxes, in integrable systems (except for 
the circle billiard) the relaxation is strongly suppressed. 
Consequently, at a given spin-orbit coupling strength, a 
chaotic cavity can be in the weak-antilocalization regime, 
whereas in a corresponding integrable cavity weak lo- 
calization will take place. Weak antilocalization is sup- 
pressed by an external magnetic field via the Zeeman 
interaction and the Aharonov-Bohm phase. The size re- 
duction of the quantum dot works against the spin re- 
laxation and the antilocalization, as well. Rashba and 
Dresselhaus interactions compensate the effects of each 
other on spin relaxation and transport. 

The degree of spin relaxation for closed orbits, which 
controls the quantum corrections to transport coeffi- 
cients, is reduced, compared to open trajectories of the 
the same length. This brings about an important differ- 
ence in the semiclassical treatment of systems with and 
without spin-orbit coupling: in the latter case, the mod- 
ulation factor due to magnetic flux can be averaged over 
arbitrary trajectories. In the limit of slow spin dynamics, 
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APPENDIX A: SPIN PROPAGATOR 

Let us consider a spin, precessing in a time-dependent 
magnetic field C(t), with the Hamiltonian H(t) = hs ■ 
C(t). The path- integral JUJl for the corresponding prop- 
agator i^(^2,Ci^) m the coherent-state representation 
yields^ 



[a*(t) - 6*(t)Ci + b(t)Q + a(t)QCi? s 



(i + ic 2 mi + icii 2 )* 



( A1 ) 

where the coefficients a(t) and b(t) are found from the 
differential equation 

^l = - l -a-C(t)W(t), W(0) = t (A2) 
for the matrix 

W ® S ( -f(t) a*% ) G SU ( 2 )- ( A3 ) 



In Eq. (|A2|) <x is the vector of Pauli matrices. As any 
SU(2) matrix, Wit) can be represented in the form: 



Wit) = 1 cos 2|) - t(T ■ m(t) sin ^ 



= e -to-ij(t)/2 



(A4) 



where m(t) = T}it)/r](t) is a unit vector. The spin s time- 
evolution operator (propagator) Kit) that belongs to the 
(2s + l)-dimensional irreducible representation of SU(2), 
satisfies the equation 



dK 
~dT 



-H(t)K = -is • C(f) K, K(0) = 1, (A5) 
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yielding 



K(t) = 



S-Jj(t) 



(A6) 



Thus, K(t), applied to the spin states at time t = 0, ro- 
tates them by the angle rj(t) about the axis m(t). The 
coherent-state propagator I|A1J) is the matrix element 

of kit)-. 

tf(C 2 ,Ci;*) = (G»l#(*)ICi) = (C 2 |e^ 8 ' T,(t) ICi)- (A7) 

Note that, setting t = (i.e., a = 1, 6 = 0) in Eq. (|A1() . 
we obtain an explicit expression for the scalar product 
(C2IC1) of two coherent states. 

It is instructive to verify Eq. (|A7J| in the case when the 
z axis is chosen along m(t). Then K{t) is diagonal in the 
|cr> basis (SecHTJ: (cr|J^|cr') = e^ 1 ^ <W . Projecting the 
spin coherent states QJ onto the basis states, 



WO = 



AS+tT 



(2s)\ 



(1 + |C| 2 ) S V (s + ay.(s-a)V 
we find from Eq. (|A7|) 

#(Ca,Ci) 

= V (CICi) a+CT (j^T (2a)! 
ff t^ s (l + lC 2 | 2 )Ml + lCi| 2 ) s (* + -a)! 

_ (i + |C 2 | 2 ) s (i + ICil 2 ) s ' 



(A8) 



(A9) 



This result agrees with Eq. (33} for a(t) = e-"'^/ 2 and 
&(*) = 0. 



APPENDIX B: LOOP CONTRIBUTION 

In Ref. the loop correction was derived for a hyper- 
bolic system with a single Lyapunov exponent A. Here we 
extend this derivation to the case with SO interaction, or, 
indeed, any interaction that does not change the classical 
trajectories in the leading semiclassical order. 

The action difference A<S(e 7i7 ') ~ Mw 2 e 2 7 ,/2A for a 
loop pair of trajectories 7 and 7' (in the absence of mag- 
netic field) depends on the crossing angle £ 7 , 7 ', as well 
as the velocity v and the mass M. Any long orbit has 
many self-crossings, each corresponding to a loop. For an 
orbit of length L, the average number of loops of length 
between I and I + dl and crossing angle between e and 
e + de is 



Pioop{l,s;L) dlde 



(L — I) sine 



dl de. 



(Bl) 



The loop length lies in the interval £ m i n (£) < I < L, 
where l m i n (s) ~ — (2v/X) ln(e/const) is the minimal pos- 
sible loop length. This length is of order of the size of 
the crossing region where the intersecting segments of 



the trajectory are interrelated by the linearized dynam- 
ics and, therefore, cannot form a closed loop. 

The loop correction to the channel-resolved transmis- 
sion coefficient |inm| 2 oop can be expressed as a sum over 
all trajectories 7(71, m) weighted with the number of 
loops l|Blll per orbit: 



loop 
pairs 



A(S(£ 75 y) 



(Mip) 



7,7' 



:2 £ \Ay 

7(n,m) 

-AS(e) 



de / <2/fl 00 p(Z,£;L 7 ) 

1' Jlmin(e) 



M<p[l - imin(e)]. 



(B2) 



By Aiip [I — Zmin(e)] we denote the average of (A / l</?) 7 , 7 ' 
over all loops with the effective length l e = I — l m i n (e). 
The contribution of the crossing region to A4(p can be 
neglected, since, in the semiclassical limit AS(e)/H ^S> 1, 
only small e are important in the integral (B2j). More- 
over, only the remaining part of the loop of length l e 
obeys the statistical assumptions under which the aver- 
age Jvltp is obtained. Substituting P\ oop (l, e; L 7 ) from 
Eq. HBlfl . we compute the I integral 



l mte (e) 



sine 



dlPi oop (l,e;L 7 )M f \y 'min (£) 



7rA c 



dl e [L 1 - l e )Mtp(l E ) 



^min 

(e) / dleMvile) 
Jo 



(B3) 



ignoring the terms of order [Z m i n (e)] 2 - When integrating 
over e in Eq. (|B2|I . we can approximate sine f=a e and 
extend the integration to infinity. Only the term propor- 
tional to l m in{e) in Eq. I|B3|) survives the integration, and 
we obtainS 



de cos 



AS(e) 



sine 
TrA r 



^min (s* ) 



[2(N + N')L e 



(B4) 



Using the results of Eqs. (B3jl and (|B4|l in Eq. (B2jl and 

applying the sum rule i|24|) , we factorize the transmission 
coefficient into a spin- and field-independent part and a 
length-averaged modulation factor (M.(p (B)) l. 



It I 2 

I L nm I 



loop 



-(N + N')- 2 LZt c / dLP L (L) / dl e M<p(l e ) 
Jo Jo 



-(N + N)~ / dl e P L (l e )Mp(l e ). 
Jo 



(B5) 



After the summation over the channels we obtain the 
loop corrections l|37|) . 
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APPENDIX C: MODULATION FACTOR FOR 
AN ARBITRARY SPIN 

To calculate the spin s modulation factor for a diago- 
nal or loop pair of trajectories, we start from Eqs. i|32|) 
and H33|) . Writing K in the exponential form (|A6|) . one 
finds 



M = Tr(K 2 ) = Tr(e 



-2is-T)\ 



E 



-1ir\a 



(Cl) 



Here, the last identity was obtained by choosing the 
z axis in the direction of r\ and using the standard matrix 
form of s z . It is sometimes convenient to represent the 
result in terms of the Gegenbauer polynomials^ Cj on 
the sphere S 3 , which generalize the Legcndre polynomials 
on S 2 . For this purpose we transform 



e~ t2va 



2s 



(-iy 



3=0 
2s 

E 

3=0 



2s-j 

E ■ 

cr=-(2s-j) 



2s _j sin [(2j + 1) g/2] 

(- 1 ■ , /0 n (G2) 

sm (77/2) 



and apply the property^ 

sin [(2j + 1)^/2] 



sin (ij/2) 

Eventually, we arrive at the expansion 



(C3) 



M = J2(-^C%(U) (C4) 
3=0 

with £ 4 = cos § [see Eq. g2J], 

In the asymptotic limit L — ► 00 of long trajectories, 
when the spin state is completely randomized, Ai (00) 
can be determined by averaging of Eq. I|C4|I over S 3 , 
which amounts to computing the scalar product of M. (£4) 
and Cq(^4) = 1 on the sphere. Then the orthogonality 
condition for the Gegenbauer polynomials^, 



(C5) 



where d£ is the surface element and 2ir 2 is the surface 
"area," yields M (00) = (— l) 2s . 



APPENDIX D: SPIN DIFFUSION 

The spin evolution in a random magnetic field is 
mapped onto a random walk on the sphere S 3 with 
initial point e 4 = (0,0,0,1). In the short-step limit 
|C| r/2 < tt/2, the probability density P d iff (£',t) to find 
the particle at point £ at time t solves the diffusion equa- 



tion 



76 



at 



777 - PV 2 P diff (£; f) = 0, P diff (£; 0) = <5(£, e 4 ), 



(Dl) 



with the diffusion coefficient D = |C| 2 r/24. Here V 2 is 
the Laplacian on S 3 , and <5(£, 64) is the 8- function on S 3 
between two points. The probability is normalized by 
Jgs d£ Pdiff (£; t) — 1, where <i£ is the surface element. 
The solution can be found in the the form of an expansion 
in the Gegenbauer polynomials^ Cj : 

-j 00 

Pm(Z;t) = Y ^J20 + 1 ) e ' j(j+2)mc H^)- (°2) 
77 3=0 

In order to determine the average modulation fac- 
tor Mdis(t), we express -M(£ 4 ) = C^iU) ~ QH&O ^ rom 
Eq. I|43|) . Taking into account the orthogonality con- 
dition l|C5|) . we compute the average for the distribu- 
tion Pdiff 

M diS (t)= f dZP m (Z;t)M(U) = 3e- 8m -l, (D3) 
Js 3 

yielding Eq. H4tj|) . This result can be generalized to an 
arbitrary spin s if one employs the representation l|C4f) : 

2.s 



The average spin polarization is, according to Eq. (|45p. 



Xdiff(<) = ^(-l) 2s - J (2j + l) e - 



(n»)di ff (*)=2K 3 3 (*)+e4 a (*)]-l- (D5) 
Using Eqs. (gSl and (gHJ), wc find 

i 4 I W = (3e-^ c l 2rt + l)/4. (D6) 

The symmetry of the problem provides 3 £ 3 2 (t) = 1 — 
£f (i), and Eq. gZJ follows. 



APPENDIX E: AREA CONTRIBUTION TO THE 
SPIN-ROTATION ANGLE 

To derive Eq. I|64|l , it is convenient to use the notation 
rj L1 = 77 in Eq. When I = 2, the BCH formula JSHJ 

applied to the product of the two segment propagators 
yields 



V 2 ,i = Vi + V2 - 2 x f)2 + C>(A X 3 Y ) 
For an arbitrary I, this result generalizes to 



(El) 



(E2) 



J=2 



where r\f\ = Y^j=iVj- The above expression can be 
proven by induction in I. Namely, we assume its validity 
for I — 1 and compute r\ l 1 from 



-IIT-T7, i/2 



(E3) 
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FIG. 18: Definition of the enclosed area [Eq. (E3)]. The 
shaded area is (1/2) Arj-i^ x Ar;. 



as in Eq. JETJ. Equation (JBiJl will follow from Eq. l(E2l) 
if the rotation angles are expressed in terms of displace- 
ments using Eqs. (|55|> and (jHSJ- The enclosed area is 
defined by 



Ae 2 = -^Ar 3 _ M x Ar,-, (E4) 
where Ar,-,i = J2Ui Ar i (Fig.[HJ>- 
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83 Equation 1191 is only valid if there are no families of direct 
trajectories (without reflections) connecting the leads—. 

84 In Ref. the diagonal contribution is defined differently 
and includes the classical part. 

85 A generalized phase-space approach for the loop contribu- 
tions, also valid for hyperbolic systems with smooth po- 
tential, is considered in Ref. [TJ 
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reversal symmetry and that C(q, p) is an odd function 
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reversal. 

87 The cyclotron radius must be larger than the system size. 

88 Semiclassical description of diffusive systems implies that 
the Fermi wavelength is much smaller than the mean free 
path. 

89 It was argued in Ref. |4^ that the saturation is due to the 
nearly periodic orbits, which have a finite measure in inte- 
grable systems. 

90 The interplay between the Rashba and the Dresselhaus 
terms in extended diffusive s yste ms was studied in Ref. [7^ 
(spin relaxation) and in Ref. |73 (weak antilocalization) . 

91 In this this case the system possesses an additional con- 
served quantity—-. 



and 



are invariant, as 



92 The partial coefficients \t nrl 
well. 

93 The requirement that the incoming and outgoing channels 
be the same was neglected. 

94 The Zeeman contribution due to spin coupling to an exter- 
nal (possibly inhomogeneous) magnetic field appears only 
in the on-site terms since it does not involve spacial deriva- 
tives. On the contrary, the SO coupling depending on the 
momentum p contributes mainly to the hopping terms. 

95 It took minutes to produce the semiclassical curves in 
Figs. 1161 and 1171 and a couple of weeks for their quantum 
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